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Abstract 

We investigate the decay of entanglement, due to decoherence, of multi-qubit systems that are 
initially prepared in highly (in some cases maximally) entangled states. We assume that during the 
decoherence processes each qubit of the system interacts with its own, independent environment. 
We determine, for systems with a small number of qubits and for various decoherence channels, 
the initial states exhibiting the most robust entanglement. We also consider a restricted version 
of this robustness-optimization problem that only involves states equivalent, under local unitary 
transformations, to the \GHZ) state. 
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I. INTRODUCTION 



Entanglement and decoherence are two intimately related phenomena that lie at the 
heart of quantum physics [l, 2, 3|. Entanglement is probably the most distinctive feature of 
quantum mechanics, its many manifestations being the focus of an intense and increasing 
research activity. From the fundamental point of view entanglement plays an important role, 
for instance, in explaining the origin of the classical macroscopic world from a quantum me- 
chanical background 3J, and also in connection with the foundations of statistical mechanics 
4j, |5j. On the other hand, the creation and manipulation of multi-partite entangled states 
have important technological applications, such as quantum computation and quantum 
metrology 6J. The phenomenon of decoherence consists, basically, of a set of effects aris- 
ing from the interaction (and concomitant entanglement-development) between quantum 

nn 

systems and their environments [2j, [3|. Almost every physical system is immersed in an envi- 
ronment and interacts with it in some way. The associated entanglement developed between 
the system and the environment leads to the suppression of typical quantum features of the 
system, such as the interference between different system's states. This constitutes the basis 
of the "decoherence program" for explaining quantum-to-classical transition 3j. 

The (internal) entanglement exhibited by a composite system undergoing decoherence 
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entanglement has 
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16[ because 



tends to decrease as the process of decoherence takes place. T 
recently attracted the interest of many researchers [3, Is,!^, 
it constitutes one of the main obstacles for the development of quantum technologies based 
on the controlled manipulation of entangled states [2J. It has been shown that in some 
cases entanglement can vanish in finite times. This phenomenon is known as entanglement 
sudden death [7, 8]. Numerous works have been devoted to the study of the robustness of 
various multipartite entangled states under the influence of different decoherence models, 
paying special attention to their scaling behavior with the size of the system jsl, [k], 11, Ijl, 
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141 ] . The dynamics of simple systems interacting with different environments has been 



experimentally studied using an all-optical device 

Some of these investigations [9| have suggested that the entanglement of multi-qubit 
systems may become more robust as the number of qubits increases. It was found that, in 
some cases, the time it takes for the entanglement of a multi-qubit system to vanish due to 
decoherence increases with the number of qubits of the system. Alas, recent results reported 
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by Aolita et al. [13j clearly indicate that, even in those cases where entanglement takes a 
long time to entirely disappear, it soon becomes too small to be of any practical use. Even 
worse, the larger the number of qubits of the system, the sooner this happens. Interestingly 
enough, the findings of Aolita et al. are fully consistent with the fast decoherence processes 
invoked (within the decoherence program) to explain the quantum-to-classic transition 3]. 

Aolita et al. [3] studied the decrease of entanglement of an A-qubits system, initially in a 
\GHZ) state, that experiences decoherence through the interaction of each of its qubits with 
an independent environment. Following this interesting and promising line of enquiry, the 
aim of the present contribution is to explore the entanglement-robustness of highly (in some 
cases maximally) entangled multiqubit-states, i.e., the ones that may be of technological 
utility. 

The paper is organized as follows. In Section II we briefly review the local decoherence 
models for multi-qubit systems that we use in the present work. In Section III we inves- 
tigate the entanglement robustness of highly entangled multi-qubit states. In Section IV 
we consider the entanglement decay corresponding to initial states equivalent under local 
unitary transformations to the \GHZ) state. In Section V we investigate whether the de- 
coherence process leads to multi-qubit states exhibiting bound entanglement. Finally, some 
conclusions are drawn in Section VI. 



II. DECOHERENCE MODELS 



The systems under consideration in the present study consist of an array on A indepen- 
dent qubits initially entangled due to a previous, arbitrary, interacting process. Each qubit 
in the composite system is coupled to its own environment; in this local environment for- 
mulation there is no communication and the entanglement between the subsystems cannot 
increase because of the locality of the involved operations. We assume that all qubits are 
affected by the same decoherence process. The dynamics of any of these qubits is governed 
by a master equation from which one can obtain a completely positive trace-preserving 
map e which describes the corresponding evolution: pi(t) = epi(0). In the Born-Markov 
approximation these maps (or channels) can be described using its Kraus representation 
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(1) 



where E~ j = 1, . . . , M are the so called Kraus operators, M being the number of operators 



needed to completely characterize the channel [17|. There are other approaches to describe 
noisy channels such as the quantum Liouville equation Qj. Using the Kraus operators 
formalism it is possible to describe the evolution of the entire N-qubit system, 



p(i) = ep(0) = '£ t Ei 1 ®...®E jN p(0) [En 



En 



(2) 



We consider the following five paradigmatic noisy channels 



1. Phase Damping 

This process describes the loss of quantum information with probability p without any 
exchange of energy. Physical examples of this process are given by the random scattering of 



a photon while traveling through a waveguide, or the pertur 
an atom when interacting with distant electrical charges 
channel are 



Dation of the electronic states in 
. The Kraus operators for this 




*„H ^_I; * = K _l (3) 



2. Depolarizing 

This one can be viewed as a process in which the initial state is mixed with a source of 
white noise with probability p. Because the channel is highly symmetric, all output states 
are unitarily equivalent. For a d-dimensional quantum system, it can be expressed as 

ep = Zl+(l-p)p (4) 

where I stands for the dxd identity matrix. The Kraus operators for this process, including 
all Pauli matrices are 

pi 



E = ^l-p'I] Et = xl^-o-i (5) 
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where p' = i = 1, 2, 3, and &i are the corresponding Pauli matrices a x , a y , a z . 
3. Bit flip, phase flip, and bit-phase flip 

These channels represent all the possible errors in quantum computation, the usual bit 
flip 0^1, the phase flip, and the combination of both, bit-phase flip. The corresponding 
pair of Kraus operators E — E 1 for each channel is given by: 



E = y/i-p/2i, E{ = VpA^; (6) 

where i = x give us the bit flip, % — z the phase flip, and i = y the bit-phase flip. Following 
Salles et.al 16|, the factor of 2 in Eq. ([6]) guarantees that at p — 1 the ignorance about the 
occurrence of an error is maximal, and as a consequence, the information about the state is 
minimum. 



III. ROBUST MAXIMALLY ENTANGLED STATES 
A. Preliminaries 

In this section we study the decay of entanglement of maximally entangled multiqubit 
pure states. There exist various measures that aim to quantify and characterize different 
features of the multipartite entanglement phenomenon. We stress that by maximally entan- 
gled states we mean those states maximizing an appropriate entanglement measure under 
the requirement that entanglement is (at least approximately) uniformly shared among all 
the system's components. These states may seem easy to characterize: all their reduced 
density matrices must be maximally mixed. A good example of a state complying with this 
requirement is the 3 qubit \GHZ) state. Alas, for systems with more than three qubits the 
\GHZ) or cat state is no longer the maximally entangled state (at least, not in the afore- 
mentioned sense). The characterization and quantification of maximally entangled states 
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N > 3 qubits has recently been the focus of an intensive research activity 
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One of the most popular measures proposed to quantify such multipartite entanglement 
is based on the use of a bipartite measure, which is averaged over all possible bipartitions 
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of the system. It is mathematically expressed by: 



[N/2] 



m=l 

^ bipart 

E (m) = j^r- £ ( 8 ) 

bipart 

Here, stands for the entanglement associated with one, single bipartition of the Af-qubit 
system. The quantity E^ gives the average entanglement between subsets of m qubits and 
the remaining N — m qubits constituting the system. The average is performed over the 
^Upart nonequivalent bipartitions, which are given by 

Kp a rt = Q if n + N/2, (9) 
NbZrt = 1(n /2 ) *n = N/2. (10) 
The total number of bipartitions is equal to 

[N/2] 

N cuts = £ Nl wart = 2 N ~ l - 1. (11) 

i=l 

Different E^ represent different entanglement features of the state. Two states may, for 
instance, share the same value of E^ and have very different values of E^ 2 \ A state may 
even maximize E^ while exhibiting very low values of E^ . It is clear that all the entangle- 
ment measures E^> have to be computed in order to describe the entanglement properties 
of a multipartite state as comprehensively as permitted by the present bi-partitions based 
approach. However, it is also desirable to try to characterize with one single number the 
"total" amount of entanglement associated with the state. Alas, any such attempt will in- 
evitably exhibit some degree of arbitrariness for the very fact that multipartite entanglement 
is a highly complex phenomenon whose description cannot be "compressed" into a single 
number. In the present work we adopted the global multiqubit entanglement measure (j7j), 
given by the average of the [N/2] different E 1 ^ associated with a state p. An alternative 
procedure would be to compute the entanglement corresponding to each of the N cuts pos- 
sible bipartitions and then evaluate directly their average, skipping the intermediate step 
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31 
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TABLE I: Number of non-equivalent bipartitions as a function of the number of qubits of the 
system N, and the dimension m of the smaller subsystem that results from the bipartition. 

of computing the quantities E^ m ' . However, this scheme implies a strong bias towards the 
most balanced bipartitions because they are more numerous than the unbalanced ones (see 
tabled]). We have performed numerical experiments using both the aforementioned ways 
of computing the global amount of entanglement of a multi-qubit state and found that the 
main conclusions reported in the present work are valid in both cases. 

We will use the negativity as our bipartite measure of entanglement because we are 
dealing with mixed states. 

Any set of unitary local operations applied to a maximally entangled state will result in a 
different state with the same amount of entanglement. When a given decoherence model is 
applied to this family of locally equivalent states quite different entanglement evolutions may 
arise. Our interest is focused on maximally entangled states that are able to retain during 
their dynamical evolution a larger amount of entanglement E than the one kept by other 
states with the same initial E— amount. Following the robustness-definition introduced in 



Ref. 24[], we will refer to these states as robust entangled ones and use the notation \^^ ob ) 
for them, where A" is the number of qubits of the system. We compute the entanglement 
dynamics of these states under the influence of the various decoherence channels introduced 
in the previous section, and compare such dynamics with that of the entanglement decay of 
the N-qubit \GHZ) and \W) states. 

To find initial states p(0) exhibiting robust entanglement under a decoherence channel 
e we use an iterative numerical search scheme akin to the well-known simulating annealing 
algorithm. During the search process a series of initial states p^(0), i — 1, 2, 3, ... of in- 



7 



creasing robustness is generated (here the word "initial" refers to the decoherence-evolution, 
not to the search algorithm). At each step of the search process, a (random) new state 
p( fc+1 )(0) slightly different from the previous one p( fc )(0) is generated. If the entanglement 
of the state ep( k+1 \0) is larger than the entanglement of ep^ k \0), the new state p( k+1 \0) 
is kept. Otherwise, the new state is rejected and a new, tentative state is generated. This 
iterative process is repeated until it converges to an initial state p(0) maximizing the entan- 
glement of the corresponding evolved state ep(0). Notice that the entanglement that we are 
maximizing is not the entanglement of p(0) itself, but the entanglement of the state ep(0) 
resulting from the action of the decoherence channel e upon p(0). Since we use the negativity 
as our measure of entanglement, the quantity to be maximized is the negativity of the mixed 
state obtained after applying the decoherence channel e to the initial state p(0). Instead of 
studying states p(0) that maximize the entanglement E[p(t)} present at a later time t, one 
could consider initial states p(0) optimizing the ratio E[p(t)}/ E[p(0)]. However, this ratio 
may adopt large values (that is, values close to one) for states with little initial (and final) 
entanglement that are not interesting as resources for quantum information processes. 



B. Present results 



As a result of the above numerical search process we were able to ascertain that some of 
the maximally entangled multiqubit states that can be found in the literature are already 
robust. However, for some systems we have found new states exhibiting a higher degree of 
entanglement robustness than the one associated with previously known states. 

At this point we begin the explanation of our novel results for the highest dimensional 
system studied in this work, that composed of 6 qubits. We find that the robust state l^c*) 
turns out to be preci sely that maximally entangled state encountered by some of the authors 



on a previous work 21] . No pure state of 6 qubits has been found that evolves to a mixed 
state with a higher amount of entanglement. In Fig. [T^i one observes that for the PD channel 
the negativity of l^ofc) ^ s almost a straight line, and a similar behavior is observed for the 
\W) state. Instead, \GHZ) quickly loses entanglement vis-a-vis the preceding examples. A 
qualitatively similar decay of entanglement is obtained in the case of the PF channel. For the 
BF-one (plotted in Fig. [lb) and also for the BPF channel, the entanglement of the \GHZ) 
is more robust than that for the preceding channels. It remains always higher than the one 
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FIG. 1: Entanglement evolution of the 6 qubits representative states under phase damping and bit 
flip decoherence models. 



associated to \W), although lower than that of |^o&)- Finally, the depolarizing channel is 
the only one for which the initial entanglement seems to be the only important factor in 
our scenario, the entanglement-decay of all our states (equally entangled initially) sharing 
such a behavior. This means that any state belonging to the set of initially maximally 
entangled states, those equivalent under LU operations to |^ ofe ), is robust according to the 
depolarizing channel. 

For lower dimensional systems the overall picture resembles that of the 6 qubits case, but 
some particularities must be pointed out. For 5 qubit systems, the robust state \^l ob ) that 
we find is not so good as the one for 6 qubits. Its entanglement decay is just maximal for 
four out of the six channels under analysis. For the two channels (BF and BPF) for which 
is not maximal its entanglement becomes lower than that of other states for large p values. 
The explicit form of l^ofe) is given by 



= -= [-|5) + |6> - |9> + |10) + 1 17) + |18) - |29> - |30) + 
V Id 

z(-|0) + |3) - |12) + |15> - |20) - |23) + |24) + |27»]. (12) 

In this equation the ket \i) denotes a member of the computational basis for a system of 5 
qubits, which can be identified by expressing the integer i in binary notation. For instance, 
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|3) is a short hand notation for |00011). 

For systems of 4 qubits, it was proved in [3] that a pure state exhibiting the theoreti- 
cally maximum amount of entanglement (that is, having all its marginal density matrices 
maximally mixed) does not exist. In Ref. [3], a promising candidate for the maximally 
entangled status was also proposed, namely, 



^ o6 ) = \HS) = -j= |1100) + |0011)+w(|1001) + |0110)) + c; 2 (jl010) + |0101)) , (13) 
with uj — — This conjecture has later received support from several numerical studies 
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23]. The entanglement decay of (fTBl resembles that of l^c*) an d its entanglement 



is always larger than that of any other state tested in our samplings. 

Contrary to what could be expected, the 3 qubit case is the most complex one. The small 
number of parameters needed to characterize a 3 qubit state is low enough to perform the 
numerical optimization process in a very short time, allowing us to identify states which 
are robust for different channels. The main difference between 3-qubit systems and the 
higher dimensional systems considered previously is that there is no 3-qubit robust state 
that simultaneously maximizes the entanglement for all the decoherence channels. All 3- 
qubit states optimizing the entanglement-robustness for a given channel quickly lose their 
entanglement when evolving under other decoherence maps. The \GHZ) provides a clear 
illustration of the peculiar features exhibited by three qubit systems. The entanglement of 
this state is robust under the action of the BF channel. However, for example, in the BPF 
channel's instance, the entanglement of the \GHZ) remains large for small p values, but 
it then decays below the expected pattern until it almost vanishes for p ~ 0.7. A similar 
behavior occurs with other channels, even when using alternative states which are robust 
for some channels but whose entanglement very quickly decays under the action of other 
decoherence processes. 



IV. DECOHERENCE OF INITIAL STATES EQUIVALENT UNDER LOCAL UNI- 
TARY OPERATIONS TO THE \GHZ) STATE 

In this Section we are going to consider the decoherence behavior of systems initially 
prepared in states equivalent (under local unitary operations) to the \GHZ) state. In par- 



10 



ticular, we are going to determine (for different decoherence channels) which one of the 
aforementioned states exhibits the most robust entanglement. 

By recourse to a numerical survey of the behavior of states equivalent to the \GHZ) one, 
we found that (among the alluded states) the one exhibiting the most robust entanglement 
for the phase damping and phase flip decoherence channels is the iV-qubits state 



\H) N = U H ® ...®U H \GHZ) 

= ^(i+r"-i-ro> (14) 

obtained by applying the Hadamard gate Uh on each of the qubits of an iV-qubits \GHZ) 
state. In equation (Tl4|) . |+) and |— ) stand, respectively, for states ^=(|0) + |1)) and ^(|0) — 

The decay of the amount of entanglement corresponding to various initial four-qubits 
states is depicted in Fig. 2 for the phase damping, depolarizing, bit flip, and bit phase 
flip decoherence channels, respectively. The four-qubits states considered in Fig. 2 are the 
\GHZ), the \H)± provided by eq. (fT4l) . as well as 1000 states generated applying random 
local unitary transformations to the \GHZ) states. We also plot the decay of entanglement 
for the robust state \HS) previously introduced in eq. (fT3|) . Excepting the \HS) state, all 
these states are equivalent under local unitary transformations to \GHZ). As we mentioned 
in the previous section, the \HS) state is the most robust 4-qubit state for all the considered 
channels. We note that for the depolarizing channel the decay of entanglement only depends 
on the initial amount of entanglement. The decay of entanglement for the \GHZ) and \H)± 
states under the PF, BF, and BPF channels is somehow equivalent, in the sense that the 
decay for the \GHZ) for the PF is the same than the decay of the \H)± for the BF. This 
is the reason why we only plot one of these channels, as the other one is equivalent. This 
behavior is easily explained because the unitary transformation Uh that maps the state 
\GHZ) into the \H)± is the same that transforms the BF channel into the PF. 

It is worth noticing that the states \H) n can be defined in terms of the recurrence relation 



\H)n 
\H)n 



_ v |o)®|fOtf-i + |i) 



\H)n-i 
\H)n-x ) 



(15) 
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FIG. 2: \H)4 vs. \GHZ) states for a 4 qubits system. In all subfigures the vertical lines correspond 
to the entanglement dynamics for 1000 states obtained from the \GHZ) by applying unitary local 
transformations. 



with 



\H) 2 = ^(|00> + |11», 
\H) 2 = ^(|01> + |10» 



(16) 



V. BOUND ENTANGLEMENT 



A multipartite entangled state has bound entanglement if it is globally entangled but is 



separable according to any of its possible bipartitions 



25] . Aolita et al. investigated the 
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FIG. 3: Entanglement evolution of different bipartitions for the robust (a) and GHZ (b) state of 
5 qubits for the depolarizing channel. The insets show a magnification of the region in which the 
entanglement vanishes. 



existence of bound entanglement in the final states resulting from decoherence processes 
acting upon initial states of the generalized \GHZ) form 13|. They found that the action 
of the depolarizing channel on a generalized \GHZ) state leads, for a certain range of p- 
values, to states for which the entanglement of the most unbalanced partitions 1 : N — 1 
vanishes while the state is still entangled according to the most balanced partitions, as for 
instance N/2 : N/2. This kind of states are said to have "bound entanglement" because no 
entanglement can be distilled from them by means of local operations, notwithstanding the 
fact that the system is globally entangled. 

We are now going to investigate the presence of bound entanglement in the iV-qubit states 
generated by the decoherence process when the initial states are the robust ones introduced in 
the previous Sections. In order to address this issue we calculate (on the states resulting from 
the action of the decoherence channel) the negativity corresponding to different partitions 
of the system. We concentrate our efforts on the depolarizing environment, because in 
this case the phenomenon of entanglement's sudden death occurs in a finite time. Fig. [3] 
shows that bound entanglement is observed both for the robust \^ r ob) and \GHZ) states. 
Actually, for all the initial states studied in this work (including a sample of 1000 randomly 
generated states) the decoherence process leads to states endowed with bound entanglement. 
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This typical behaviour indicates a remarkable capacity of the environment to create bound 
entanglement. However, special initial states exhibit interesting differences when the decay 
of entanglement associated with different types of bipartitions is considered. For the robust 
states \^rob) introduced in this work, as well as for all the random states that we have 
generated, the entanglement of the most balanced partitions N/2 : N/2 vanishes sooner 
than the entanglement corresponding to the most unbalanced partitions 1 : N — 1. This 
is in contraposition to what happens with the generalized \GHZ) states, as was found by 
Aolita et al. The \ W) and the \H) n states also share this uncommon behaviour: the p- values 
for which the entanglement of the most unbalanced partitions 1 : N — 1 vanish corresponds 
to states that are still entangled according to more balanced partitions. 



VI. CONCLUSIONS 



In the present effort we investigated the decay of the amount of entanglement of a multi- 
qubit system experiencing a decoherence process. We considered models of decoherence 
characterized by independent environments respectively interacting with each of the system's 
qubits. We performed a systematic numerical search of the initial pure states exhibiting the 
most robust entanglement under decoherence scenarios described by five different channels: 
phase damping, depolarizing, bit flip, phase flip, and bit-phase flip. For systems of 4 qubits 
we found that the state having the most robust entanglement is the \HS) state introduced 
by Higuchi and Sudbery 18|, which was conjectured to correspond to a global maximum 
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221 ] . In the case of six qubits the 
ready known state of 
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221 ] . The behavior 



of entanglement for four-qubits systems 

state maximizing entanglement robustness again coincides with an a 
maximum entanglement that had been the subject of previous studies 
of the states optimizing robustness for 4 and 6 qubits is similar: both states maximize 
entanglement-robustness for each of the six decoherence channels considered in the present 
work. 

We also were able to determine a five-qubits state with highly robust entanglement. 
However, this state exhibits optimal robustness only for four of the decoherence channels 
considered here. Under the effect of the remaining decoherence channels, and for large 
enough values of the "time" parameter p, this initial state evolves into a state with less 



entanglement than that obtained from other initial states. When this happens, however, the 
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entanglement exhibited by the system (regardless of the initial state) is already too small 
to be of any practical use. 

Unexpectedly, the case of three-qubits turned out to be the most complex one. Our 
numerical results indicate that for three-qubits systems there is no state simultaneously 
maximizing the entanglement-robustness for all (or, at least, for most of) the decoherence 
channels studied here. In other words, each state exhibiting optimal robustness for one of 
the channels was found to be non optimal for some of the other channels. 

We also discussed a restricted version of the robustness' optimization problem involving 
only states equivalent to the \GHZ) one under local unitary transformations. A state \H)n 
obtained by applying the Hadamard gate to each of the qubits of a multi-qubit system in 
the \GHZ) state was found to play an interesting role. This special state maximizes the 
entanglement-robustness for four decoherence channels. 

Finally, we studied the phenomenon of bound entanglement determining whether initially 
entangled multi-qubit states evolve (as the decoherence process takes place) into states 
characterized by bound entanglement. We focused our attention on the depolarizing channel, 
this being the channel leading to entanglement sudden death in a finite time. We found that 
the ability of the environment to create bound entanglement is a quite general feature that 
appears for all the states studied in this work. 
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